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Abstract. Based on 14 data points obtained with near-infrared speckle interferometry and covering an almost entire revolution, 
we present a first visual orbit for the low-mass binary system Gliese 22 AC The quality of the orbit is largely improved with 
respect to previous astrometric solutions. The dynamical system mass is 0.592 ± 0.065 M@, where the largest part of the error is 
due to the Hipparcos parallax. A comparison of this dynamical mass with mass-luminosity relations on the lower main sequence 
and theoretical evolutionary models for low-mass objects shows that both probably underestimate the masses of M dwarfs. A 
mass estimate for the companion Gliese 22 C indicates that this object is a very low-mass star with a mass close to the hydrogen 
burning mass limit. 

Key words. Stars: individual: Gliese 22 — Stars: binaries: visual — Stars: low-mass, brown dwarfs — Techniques: high 
angular resolution 



1. Introduction 

M dwarfs are the dominant population of the Galaxy as well 
in numbers as in stellar mass contribution. Furthermore, they 
mark the transition regime between stars and the now well es- 
tablished classes of substellar objects. Given these important 
properties, it is problematic that there is only a small number 
of empirically determined stellar masses at the lower end of the 
main sequence. In addition, most of these dynamical masses 
are affected with large uncertainties (e.g. Henry 1998). This 
means that the mass-luminosity relation for M dwarfs is not 
well calibrated. Moreover, theoretical evolutionary models for 
low-mass objects that cover also the substellar regime are not 
well checked with dynamical masses. 

As a contribution to a solution of these problems we are car- 
rying out a program aiming at a determination of visual or- 
bits and thus dynamical masses for M dwarf binaries. Speckle 
interferometry with array cameras in the near infrared allows 
highly precise measurements of the relative astrometry in sub- 
arcsecond binary or triple systems. 

This program has already led to orbit determinations for the 
very low-mass systems Gliese 866 (Woitas et al. 2000) and 
LHS 1070 BC (Leinert et al. 2001). The companions in these 
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systems appear to have masses close to the stellar/substellar 
limit at M « 0.075 M . In this paper we discuss a first visual or- 
bit determination for Gliese 22 AC that was briefly announced 
in IAU Commission 26 Information Circular 147 (Docobo et 
al. 2002). 

A visual companion to the M2 star Gliese 22 (other desig- 
nations: HIP 2552, BD +66°34, ADS 440, V 547 Cas) was 
first reported by Espin & Milburn (1926). At the time of this 
detection the projected separation was 2'.'79. The orbital mo- 
tion of this companion with a period of as 320 yr has been 
monitored since its detection, and most recent orbital elements 
are given by Lampens & Strigachev (2001). Alden (1947) 
found that the primary component of this pair is itself an as- 
trometric binary. Herafter, we will refer to this close pair as 
Gliese 22 AC and to the more distant third component as 
Gliese 22 B. Hershey (1973) presented orbital elements for 
Gliese 22 AC based on a rich collection of astrometric plates 
from the Sproul Observatory. This calculation was refined by 
Heintz (1993) adding more data points and Soderhjelm (1999) 
including the Hipparcos parallax of Gliese 22 into the analy- 
sis. The first two resolved visual observations of Gliese 22 AC 
were reported by McCarthy et al. (1991) using near-infrared 
speckle interferometry. Since the epoch of these first two mea- 
surements we regularly observed this pair obtaining 12 more 
data points, which uniformly cover almost an entire revolution. 
Based on these data, we present in this paper a visual orbit and 
a dynamical system mass that are more precise than the previ- 
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ous astrometric solutions. We will describe the techniques of 
observations and data analysis in Sect. 2 and present the result 
of the orbit calculation in Sect. 3. In Sect. 4 we will discuss im- 
plications of the derived dynamical system mass on the mass- 
luminosity relation and theoretical models for very low-mass 
objects. 

2. Observations and Data Analysis 

The database for our visual orbit determination for 
Gliese 22 AC is given in Table 1 . The observations num- 
bered 1 and 2 were taken from McCarthy et al. (1991) while 
the other measurements are published here for the first time. 
Observations 3 and 4 made use of one-dimensional speckle 
interferometry. This observing technique and the reduction of 
these data are described in Leinert & Haas (1989). 
All other data points have been obtained with the near-infrared 
cameras MAGIC and OMEGA Cass at the 3.5-m telescope 
on Calar Alto. Both instruments are capable of taking fast 
sequences of short time exposures (f exp ~ 0.1 s) and in this 
way allow speckle interferometry with two-dimensional 
detector arrays. Typically we have taken 1000 short expo- 
sures for Gliese 22 and the nearby PSF calibrator (single star) 
SAO 1 1358 in the K band (A = 2.2 ^m). A detailed overview of 
the data reduction and analysis has been given by Kohler et al. 
(2000). Briefly, we obtain the modulus of the complex visibil- 
ity by deconvolving the power spectrum of Gliese 22 with that 
of the PSF calibrator. The phase is recursively reconstructed 
using the algorithm by Knox & Thompson (1974) and the 
bispectrum method (Lohmann et al. 1983). As an example we 
show in Fig. 1 modulus and bispectrum phase obtained from 
the observation at 3 Nov 2001. In the two-dimensional obser- 
vations Gliese 22 B is also in the detector array. Therefore, 
we fit a triple star model to the complex visibility and in this 
way determine the relative astrometry and the flux ratios of 
the companions B and C with respect to Gliese 22 A. Pixel 
scale and detector orientation are derived from astrometric fits 
to images of the Orion Trapezium cluster core where precise 
astrometry has been given by McCaughrean & Stauffer (1994). 
These calibration observations are however only available for 
the observations since 1995 (No. 9 to 14 in Table 1). For the 
previous observations we have determined pixel scale and 
detector orientation from visual binary stars that either have 
well known orbits (a Psc, £ Aqr) or show no measurable 
orbital motion within 10 yr (RNO 1 BC). Position angles and 
projected separations of these binaries were calibrated with 
the help of the Trapezium cluster in subsequent observing 
runs. As an example we show our Trapezium-calibrated 
measurements of a Psc in Table 3, together with the prediction 
of the ephemerides from Scardia (1983). In this way all our 
two-dimensional speckle observations were placed into a 
consistent system of pixel scale and detector orientation. 
Position angles, projected separations and flux ratios for the 
companions C and B are given in Tables 1 and 2. The formal 
uncertainties of the relative astrometry are typically 5-10 
milli-arcsec in x and y for the 2D data. These errors seem to 
be reasonable since they are in the same order of magnitude as 
the residuals of the orbital fit (see Table 3). 



Table 5. Ephemerides for Gliese 22 AC 



Epoch 


[deg] 


P["] 


2003.0 


351.9 


0.485 


2004.0 


8.5 


0.521 


2005.0 


24.2 


0.517 


2006.0 


41.1 


0.489 


2007.0 


60.3 


0.455 


2008.0 


82.1 


0.435 


2009.0 


104.7 


0.441 


2010.0 


125.5 


0.469 


2011.0 


143.4 


0.506 


2012.0 


159.2 


0.531 


2013.0 


174.3 


0.527 


2014.0 


190.8 


0.483 


2015.0 


212.7 


0.400 



Table 6. (O-C) residuals for our visual orbit for Gliese 22 AC and the 
most recent astrometric orbit for this system presented by Soderhjelm 
(1999). 6 is given in degrees and p in arcseconds. 



Epoch 


Soderhjelm (1999) 


This work 




6 


P 


e 


P 


1989.778 


+2.4 


+0.006 


+3.5 


-0.041 


1989.939 


+0.1 


+0.010 


+ 1.8 


-0.033 


1990.693 


-3.3 


-0.026 


+0.6 


-0.055 


1991.717 


-13.2 


+0.032 


-6.8 


+0.007 


1993.745 


-14.7 


+0.056 


-0.3 


0.000 


1994.074 


-16.2 


+0.065 


+0.1 


+0.001 


1994.704 


-19.5 


+0.085 


+0.6 


+0.010 


1994.947 


-21.1 


+0.102 


+0.7 


+0.024 


1995.767 


-26.9 


+0.082 


0.0 


-0.003 


1996.739 


-32.7 


+0.083 


0.0 


+0.003 


1997.876 


-38.7 


+0.020 


-0.6 


-0.023 


2001.107 


-6.3 


-0.070 


+0.2 


+0.004 


2001.838 


-1.6 


-0.012 


-0.8 


+0.009 


2002.797 


-1.5 


-0.009 


0.0 


-0.010 



3. Results 

To calculate the visual orbit indicated in Fig. 2, we used the 
method proposed by Docobo (1985). A weight from 4 to 10 
was assigned to the individual measurements (see Table 1) to 
take into account seeing conditions and the quality of different 
instruments. The resulting orbital elements and their uncertain- 
ties are presented in Table 4 together with the results of the 
three previous astrometric orbital solutions. Ephemerides for 9 
and p until epoch 2015.0 are given in Table 5. The residuals 
of our data points with respect to our visual orbit and the most 
recent astrometric orbit by Soderhjelm 1999 (Table 6) indicate 
that our visual orbit represents a significant improvement. This 
is also evident from Fig. 2 (right panel) where our data points 
are plotted together with all four orbits. The precision of the or- 
bital elements is rather high, and we expect only minor changes 
to this orbit. Thus, grade 2 (good orbit) can be assigned to it ac- 
cording to the grading scheme described in the Sixth Catalog 
of Orbits of Visual Binary Stars (Hartkopf & Mason 2003). 
The system mass derived from a and P, using the Hipparcos 
parallax (98.74 ± 3.37 milli-arcsec) is equal to 0.592 ± 
O.O65M . This represents a rather high relative accuracy as the 
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Fig. 1. Modulus (left) and bispectrum- 
phase (right) for Gliese 22, derived from 
1000 short (r = 0.2 s) exposures obtained at 
A = 2.2 with OMEGA Cass at the 3.5- 
m telescope on Calar Alto on 3 November 
2001. North is to the top and east to the 
left. The broad stripe pattern in the mod- 
ulus is due to Gliese 22 AC while the fine 
patterns (roughly from top to bottom) are 
produced by the more distant companion 
Gliese 22 B. The phase has been cut off at 
the Nyquist frequency that is 5.3 arcsec" 1 
for the adopted pixel scale of 07095/pixel. 

Table 1. Overview of all observations of Gliese 22 AC used for the visual orbit determination in this paper. Except for the first two data points 
all observations have been carried out at the 3.5-m telescope on Calar Alto. 



No. 


Date 


Epoch 


Position angle 


Projected 


Weight 


Filter 


Flux ratio 


Instrument 








angle [°] 


separation [mas] 






FclF A 


(or reference) 


1 


12.10.1989 


1989.7803 


42.7 ± 2.5 


451 ± 20 


4 


K 


0.167 ±0.009 


McCarthy et al. 1991 


2 


10.12.1989 


1989.9418 


43.9 ± 2.5 


453 ± 20 


4 


H 


0.143 ±0.008 


McCarthy et al. 1991 


3 


11.09.1990 


1990.6954 


57.1 ± 6.3 


405 ± 39 


5 


K 




ID 


4 


20.09.1991 


1991.7201 


71.6+ 1.6 


443 ± 13 


5 


K 




ID 


5 


30.09.1993 


1993.7474 


122.5 ± 0.5 


464 ± 11 


7 


K 


0.151 ±0.004 


MAGIC 


6 


28.01.1994 


1994.0767 


129.3 ±4.3 


478 ± 40 


10 


K 




MAGIC 


7 


15.09.1994 


1994.7064 


141.1 ±0.4 


510 ±5 


9 


K 


0.141 ±0.007 


MAGIC 


8 


13.12.1994 


1994.9501 


145.2 ± 0.2 


532 ± 5 


9 


K 




MAGIC 


9 


08.10.1995 


1995.7693 


157.5 ± 1.7 


526 ± 16 


9 


K 


0.158 ±0.004 


MAGIC 


10 


28.09.1996 


1996.7448 


172.1 ±0.2 


533 ±4 


10 


K 


0.141 ±0.004 


MAGIC 


11 


17.11.1997 


1997.8789 


190.2 ± 0.3 


460 ±4 


7 


K 


0.171 ±0.008 


MAGIC 


12 


09.02.2001 


2001.1095 


302.3 ± 0.3 


334 ±4 


9 


K 


0.185 ±0.005 


OMEGA Cass 


13 


03.11.2001 


2001.8406 


325.5 ± 0.5 


402 ± 6 


8 


K 


0.151 ±0.011 


OMEGA Cass 


14 


19.10.2002 


2002.7994 


348.2 ± 0.4 


463 ±5 


9 


K 


0.150 ±0.005 


OMEGA Cass 




Table 2. New observations of the wide pair Gliese 22 AB obtained at the 3.5-m telescope on Calar Alto. 



No. 


Date 


Epoch 


Position angle 


Projected 


Filter 


Flux ratio 


Instrument 








angle [°] 


separation 




FbIF a 




1 


30.09.1993 


1993.7474 


164.4 ± 0.5 


4'.'24 ± 0T09 


K 


0.199 ±0.004 


MAGIC 


2 


28.01.1994 


1994.0767 


165.3 ± 0.6 


4'.'28 ± 0706 


K 




MAGIC 


3 


15.09.1994 


1994.7064 


167.1 ± 0.4 


4'.'29 ± 0704 


K 




MAGIC 


4 


13.12.1994 


1994.9501 


166.9 ±0.1 


4'.'29 ± a/oi 


K 


0.267 ± 0.004 


MAGIC 


5 


08.10.1995 


1995.7693 


167.4 ± 0.2 


4729 ± 0701 


K 




MAGIC 


6 


28.09.1996 


1996.7448 


168.7 ± 0.1 


4728 ± 0701 


K 


0.245 ± 0.004 


MAGIC 


7 


17.11.1997 


1997.8789 


170.0 ± 0.3 


4723 ± a/01 


K 


0.325 ± 0.004 


MAGIC 


8 


09.02.2001 


2001.8406 


172.8 ± 0.2 


4700 ± 0701 


K 


0.246 ± 0.006 


OMEGA Cass 


9 


03.11.2001 


2001.8406 


172.8 ±0.1 


3799 ± a/01 


K 


0.270 ± 0.009 


OMEGA Cass 


10 


19.10.2002 


2002.7994 


173.0 ±0.1 


37956 ± a/005 


K 


0.298 ± 0.005 


OMEGA Cass 



direct dynamical mass sum determination concerned. While a 
good quality orbit (in particular, a and P are determined with 
an accuracy of almost 1%) is obtained, the accuracy of mass 
determination is still 11%. The latter is a good accuracy itself 
since the usual values for post-Hipparcos epoch are in the order 
of 20 - 25% (see Martin et al. 1998, Soderhjelm 1999) but it is 
worth noting that it could be even better. 
The principal reason in this (and many other) cases is the low 
accuracy of the Hipparcos parallax (more than 3%) whose con- 



tribution to the mass error is 88% while semimajor axis and 
period contribution are almost insignificant (5% and 7% re- 
spectively). Gliese 22 AC is thus a good example of the over- 
all mass accuracy deterioration due basically to the low rela- 
tive accuracy of the parallax. Therefore, one must state that the 
sensibly higher accuracy next generation post-Hipparcos paral- 
laxes are needed to drastically improve the relative accuracy of 
direct mass determination. Should semimajor axis, period and 
parallax be each determined with 1 % accuracy, a mass deter- 
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Table 3. Our observations of a Psc that employ an astrometric calibration derived from images of the Trapezium cluster core, compared to 
the prediction of the ephemerides from Scardia (1983). The mean residuals were used as a correction for the ephemerides. In this way the 
measurements of Gliese 22 before 1995 were put into a consistent system of pixel scale and detector orientation. 



Epoch 



Instrument Telescope 



Measured 
(calibrated with Trapezium) 
PA [deg] d [arcsec] 



Predicted 
(Scardia 1983) 
PA [deg] d [arcsec] 



APA 



Residuals 



Ad 



1995.526 ESO NTT 

1996.641 ESO 3.6-m 

1997.8843 CA 3.5-m 

1999.6681 CA 3.5-m 



SHARP II 275.8 ± 0.1 1.848 ± 0.005 276.5 1.868 

ADONIS 274.5 ± 0.1 1.870 ± 0.008 275.8 1.861 

MAGIC 274.8 ± 0.2 1.838 ± 0.012 275.0 1.853 

OMEGA Cass 273.1 ± 0.1 1.852 ± 0.004 273.8 1.842 



-0.7 
-1.3 
-0.2 
-0.7 



-0.020 
0.009 
-0.015 
0.010 



mean residuals -0.7 ± 0.2 -0.004 ± 0.008 



Table 4. The orbits for Gliese 22 AC 



Element 


Hershey (1973) 


Heintz (1993) 


S6derhjelm(1999) 


This work 


P(yr) 


15.95 ± 0.22 


16.0 


15.4 


16.12 ±0.2 


T 


1956.0 ± 2.8 


1988.9 


1989.0 


2000.47 ± 0.20 


e 


0.05 ± 0.07 


0.0 


0.0 


0.18 ±0.03 


al"] 


0.51 


0.525 


0.450 


0.529 ± 0.005 


i [deg] 


45.0 


42.0 


27.0 


46 ± 1 


Q. 


167.0 


171.0 


24.0 


179.7 ± 1 


w 


160.0 


0.0 


0.0 


93.0 ± 5 




Fig. 2. Upper panel: Our visual orbit for Gliese 22 AC, together with the residuals of the individual data points. The line of nodes is also 
indicated. The motion of the companion is anticlockwise. Lower panel: Comparison of our data points with all four orbits from Table 4: Our 
orbit (thick solid line) and the astrometric orbits from Hershey (1973, dot-dot-dashed), Heintz (1993, dashed) and Soderhjelm (1999, dotted). 



mination accuracy of 5% can be achieved which is especially 
important for the lower end of luminosities in the HR diagram. 



4. Discussion 

We can now compare our dynamical mass with predictions 
from the present mass-luminosity relations on the lower main 
sequence as well as from theoretical evolutionary models for 
low-mass objcts. The absolute (system) K band magnitude for 
Gliese 22 AC is 6.26 + 0.12 mag (McCarthy et al. 1991) The 



Woitas et al.: Visual orbit for Gliese 22 AC 



5 



mean flux ratio from all measurements in Table 1 is Fc/Fa = 
0. 156 ± 0.005. This results in absolute magnitudes for the com- 
ponents: 

K A = 6.42 ±0.12 

K c = 8.43 ±0.13 (1) 

The K band mass-luminosity relation given by Henry & 
McCarthy (1993, their Eq. 2) translates this into component 

masses: 

M A = O.378!°;°gM 

M c = 0.1 36 + _° Z M O (2) 

The sum of these masses is 0.514 ± 0.029 Mq. The dy- 
namical system mass of 0.592 ± 0.065 Mq is in line with 
this result only within its 2cr error. McCarthy et al. (1991) 
did a similar calculation and found that the mass sum de- 
rived from the infrared magnitude-mass relations of Henry & 
McCarthy (1990) and the dynamical mass were remarkably 
consistent. For the latter they used however the orbit calcu- 
lation by Hershey (1973) which is worse than our visual or- 
bit (see Fig. 2) and results in a probably too low system mass 
of 0.485 ± 0.071 M Q . Our result indicates that the the K band 
mass-luminosity relations by Henry & McCarthy probably un- 
derestimate stellar masses for M dwarfs. A similar result can 
be obtained for the optical mass-luminosity relation given by 
Henry et al. (1999). Using resolved photometry taken with the 
HST Fine Guidance Sensors and updated parallax information 
(van Altena et al. 1995), these authors derive a system mass of 
0.489 ± 0.041 M G . This is again distinctly lower than our dy- 
namical value. 

The components' K band magnitudes from Eq. 1 can also be 
used to derive masses from the theoretical evolutionary models 
for low-mass stars and substellar objects presented by Baraffe 
et al. (1998, 2002). For this purpose one has to assume an 
age for Gliese 22. Since this system shows no signs of chro- 
mospheric activity (Ha is in absorption, see Herbst & Miller 
1989), we adopt a lower age limit of 10 8 yr. As can be seen 
from Table 7, the system mass derived from the theoretical 
models is < 0.54 Mq for all ages above this value. This is again 
less than the dynamical system mass of 0.592 Mq, but compa- 
rable to this empirical result within the uncertainties. The re- 
striction of this discussion to ages > 10 8 yr causes no bias since 
lower ages would yield much lower (and thus unrealistic) mass 
estimates. 

Gliese 22 AC is a variable star that shows flare events with am- 
plitudes of 0.6 mag at optical wavelenghts (Pettersen 1975). 
One may ask how this property influences the previous dis- 
cussion. Even in the unlikely case that the components' mag- 
nitudes from Table 1 were affected by flares, our conclusions 
would not be altered. Temporary higher luminosities would 
lead to higher mass estimates, and thus the masses inferred 
from the mass-luminosity relation and theoretical evolutionary 
models would be even lower than our values if derived from 
observations in the quiescent state of a flare star. 



Table 7. Masses for Gliese 22 A and C derived from their K magni- 
tudes (Eq. 1) using the theoretical evolutionary models from Baraffe 
etal. (1998, 2002). 

Age[yr] M A [M Q ] M C [M ] M S> „[M ] 



10 s 


0.36 


0.095 


0.455 


10 85 


0.40 


0.14 


0.54 


10" 


0.39 


0.14 


0.53 


10 10 


0.39 


0.14 


0.53 



The lower mass component Gl 22 C is apparently a very 
low-mass star with a mass close to the hydrogen burning mass 
limit atM s 0.075 M . Therefore, it would be very interest- 
ing to derive its individual dynamical mass from an absolute 
orbit and compare it to the predictions of theoretical models. 
It will probably indeed be possible to disentangle the relative 
orbit for Gliese 22 AC into two absolute orbits for both compo- 
nents using the distant companion Gliese 22 B as astrometric 
reference. However, to obtain most reliable absolute orbits, a 
complete coverage of an orbital revolution of the AC pair is 
strongly desirable, and this is thus beyond the scope of this pa- 
per. 

5. Summary 

With the first calculation of a purely visual orbit from inter- 
ferometric data we have largely improved on the accuracy of 
the orbital elements for the low-mass binary Gliese 22 AC and 
derived a dynamical system mass of 0.592 ± 0.065 Mq. The 
uncertainty of 1 1% is mostly due to the error of the Hipparcos 
parallax, while future (p, 8) measurements will probably result 
in only minor changes of the orbital parameters. Based on the 
K band magnitudes of the components we have estimated their 
masses also from the mass-luminosity relation given by Henry 
& McCarthy (1993) and from the theoretical evolutionary mod- 
els for low-mass objects by Baraffe et al. (1998, 2002). In both 
cases the obtained mass sum is lower than the dynamical sys- 
tem mass. The component Gliese 22 C is apparently a very 
low-mass star with a mass around 0. 1 Mq and is thus located 
at the very end of the lower main sequence. 
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